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Abstract 



Finite-dimensional reductions of the 2D dispersionless Toda hierarchy constrained by the 
"string equation" are studied. These include solutions determined by polynomial, rational 
or logarithmic functions, which are of interest in relation to the "Laplacian growth" or Hele- 
, Shaw problem governing interface dynamics. The consistency of such reductions is proved, 

J> ■ and the Hamiltonian structure of the reduced dynamics is derived. The Poisson structure of 

OO I the rationally reduced dispersionless Toda hierarchies is also derived. 

o 

^ ; 1 Introduction. 



I This paper concerns rational and logarithmic reductions of the 2D dispersionless Toda hierarchy of inte- 

. grable equations (henceforth 2dToda). The subject is motivated by important applications to problems 

' in interface dynamics and statistical physics. 

Laplacian growth is a process that governs the dynamics of the boundary in the plane separating two 
i^H ' disjoint, open regions 2?+ and X'_ in which harmonic (scalar) fields are defined. These may be interpreted 

. as the pressure fields for two incompressible viscous fluids (Hele-Shaw problem). The movement of the 

boundary is determined (according to Darcy's law, in the case of viscous fluids) by equating the normal 
velocity of the boundary to the boundary value of the gradient of the field. In particular, one region (say, 
the "interior" region 2?+ ) may be chosen to be bounded and have constant harmonic field (corresponding 
to zero viscosity) with the boundary condition for the "exterior" I?_ region at infinity to be such that 
$H ' there is a unit sink, implying that the area of the interior region grows linearly in time [16]. Denoting 

the harmonic field (e.g. the pressure) in the exterior region by P{X,Y), this satisfies the conditions: 



AP{X,Y) = 

P ^ (47r)-iln(X2 + y2) as X^ + Y^^oo 



2 ,v2 ^ ^ (1) 



in the cartesian coordinates (X,Y). The normalized exterior normal velocity at the boundary is given 
by: 

= -nVP. (2) 

where n denotes an outward normal to dT>^ In the case where the boundary is an analytic curve it 
is usual to use the Riemann mapping theorem to introduce a time-dependent conformal map from the 
exterior of the unit circle in the complex w plane to the exterior region P_ in the "physical" plane 
z = X + iY taking the unit circle to the boundary 

z = z{w,x), w = exp(z0), < < 27r, (3) 

where x stands for the physical time of the Hele-Shaw problem. We choose this unusual notation for 
time for consistency with that used in the literature on the dispersionless integrable systems. 
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Simple considerations [15], [4], [12] show (sec Appendix A) that equations (1, 2) are equivalent to the 
following equation 

J / dz dz\ f dz{w,x) dz{l/w,x) dz{w,x) dz{l/w,x)\ ^ 
\d(j)dx ) \ dw dx dx dw J ' 

where bar stands for complex conjugation (and w = on the boundary curve). In our notations 
^i'^) = J2i ^i'^'' if -^M = J2i while z{w) = Y,i ZiW\ 

Known as the Galin-Polubarinova equation in the Hele-Shaw problem, eq (4) plays an essential role in 
the theory of infinite-dimensional integrable hierarchies in the dispersionlcss limit. The relation between 
the boundary dynamics above and the dispersionlcss limit of the integrable Toda hierarchy constrained 
by (4) was shown in [12]. 

Equation (4) may be interpreted as a constraint on an infinite commiiting set of dynamical systems 
defined in the space z{x^w) of one-parameter families of conformal maps. This constraint represents 
fixed points of an "additional symmetry" [13] and is called the "string equation" in the theory of in- 
tegrable systems. The most interesting aspect of such constrained 2dToda flows is that they admit 
finite-dimensional reductions, which include so-called "multi-finger" solutions [7]. These solutions are of 
great importance in practical applications and describe numerous phenomena, such as viscous fingering 
in a Hele-Shaw cell [7], [10], [17] and pattern formation in the quantum Hall effect [1]. 

In what follows, we consider finite-dimensional reductions of (4) in the context of the 2dToda hier- 
archy. We first study algebraic solutions of the problem, ignoring the real structure and treating 2;, ^ as 
independent functions, and w as a formal variable. Returning to the applications to interface dynamics 
we identify bar with complex conjugation. 



2 2dToda hierarchy and string equation. 

The 2dToda hierarchy is a dispersionlcss limit [18] of the two-dimensional Toda hierarchy and is defined 
in terms of two functions z{'w,x) and z{'w~^,x) of the form: 



z{w,x) = r{x)w + '^^Uk{x)w , (5) 

fe=0 

00 

z{w~^ , x) = r{x)w~^ +'^^^Uk{x)w'^ . (6) 



fe=o 



The 2dToda flow equations are 



dtkZ = {Hk, z] dt^z = {Hk, z} , . 

dt,J = {Hk,z} dt^z = {Hk,z} 



where the " Poisson-Lax" bracket notation here denotes 

{/) 9} ^ ^ g^ (8) 

and the coefficients r{x),Uk{x),Uk{x) are viewed as coordinate functions on the phase space. The evolu- 
tion functions are defined as follows 

Hk = {z% + l/2{z'')o, Hk = iz'')- + l/2{z%, (9) 

where subscripts ±, denote the negative/positive and zero parts of the formal Laurent expansion in w 

It is important to note that, despite the notation, (8) is not really a Poisson bracket defined on 
the infinite-dimensional phase space of the 2dToda system (7) (with coordinates r,Uk,Uk), but rather a 
"quasiclassical" — > limit of commutators in the Lax representation of the dispersive version of (7) 
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(where h stands for the lattice spacing, see [18]). Values (9) are dispersionless analogs of the upper or 
lower diagonal parts of the powers of Lax matrices. 

The 2dToda system is nevertheless an integrable Hamiltonian system of PDEs for the functions 
r{x,ti, ...,ti...), Uk{x,ti, ...,ti...), Uk{x,ti, ...,ti...), obtained by equating coefficients of (7) as Laurent 
polynomials in the dummy variable w. One can easily verify that vector fields (7) commute, i.e 

( z\_ ( z\ z\^_d^( z\ ( f M 

dtidtj \ z J dtjdti \z)' diidtj \z) dtjdU \ z J ' dtjdU \ z J diidtj \ z J 

due to the "zero curvature" conditions 

9Hi dH, dH, dHj - dSi dHj - - 

which follows from equations of motion (7) and definitions (8), (9) . 

The Galin-Polubarinova equation (4) for the Hele-Shaw problem written in terms of the Poisson-Lax 
bracket (8) 

{z{w,x),z{w~'^,x)} = 1 (10) 

is the string equation. 

It is fundamental that equation (10) is invariant under the flows generated by (7). Indeed 

dt, {z, z} = {dt,z, z) + {z, dt^z) = {{Hk, z}, z} + {z, {Hk, z}} = {{z, z}, Hk} = 0, (11) 

where we have used (7), (10) and the Jacobi identity for the Lax-Poisson bracket (8). 

Thus the string equation (10) defines an invariant manifold under the Toda flows (7) and a reduction 
of the Toda hierarchy. On the other hand, the Toda flows may be viewed as symmetries generating new 
solutions of Eq (10). 

3 Reductions of the 2dToda hierarchy constrained by the string 
equation. 

The reduction of the 2dToda hierarchy by the string equation is still a compatible set of infinite- 
dimensional dynamical systems. 

Indeed, as seen from (5)-(6) the string equation (10) is a system of ODE's of the form 

dr / dx — R{r,Ui, ...,ui, ...), dui/dx = Ui{r,ui, ...,ui, ...), dui/dx = Ui{r,ui, ...,ui, ...) (12) 

In what follows we will be interested in further "functional" reductions where z, z are polynomial, 
rational or logarithmic functions of w. As shown below, such reductions are consistent with (7) (i.e. the 
corresponding ansatz for z is preserved by the 2D Toda flows) if the string equation (10) holds. Thus, 
for consistency we need a double ("functional" plus "string") reduction. These pairs of reductions define 
finite-dimensional invariant sub-manifolds in the phase space of the general 2dToda hierarchy. Indeed, 
functional reduction leaves a finite number of discrete indices in the ansatz for z as a function of w, 
while the string equation fixes the dependence of z on x, leaving finite number of degrees of freedom. 
These degrees of freedom are connected with the integration constants of system (12), which becomes 
finite-dimensional after functional reductions. 

3.1 Polynomial reductions. 

We begin with polynomial reductions of the 2d Toda hierarchy 



z{w) = rw + Ujw' (13) 

i=-N 
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N 

z{'w~^) = rw~^ +'^Uiw'- (14) 

i=0 

constrained by the string equation (10). 

The following proposition states the consistency of the polynomial reductions under the Toda flows: 

Proposition 1 If the string equation (10) holds, then (13), (14) is a finite dimensional subspace (of the 
phase space (5), (6)) invariant under the ti and ii 2dToda flows, < i < N + 2 (7). This subspace has 
dimension 2N + 3 with local coordinates chosen as the initial values of the solutions r = r{x), Ui = Ui{x), 
Ui = Ui{x) to the reduced string equation (12). 

Proof. We must prove that z remains of the form (13), under the flows generated by H^, Rk provided 
the string equation (10) holds. In other words, we have to show that the evolution does not change the 
highest and lowest degrees of the Laurent polynomial (13). The proof for z is analogous. 

1. First we proceed with evolution of z generated by H^. The lowest degree of w in ^; (13) is —N. 
Since 

Hk = (z'^)>o + l/2(z'^)o = hk{x)w'' + hu-i{x)w''-^ + ... + ho{x) (15) 

is a Laurent polynomial of positive degree in w, the lowest degree of the bracket {z, Hk} is not less than 
—N, as follows from (8). On the other hand, the complement z'^ — Hk of (15) is a polynomial in l/w 
and so, 

{Hk,z} = -{z''~Hk,z} 

is a Laurent polynomial whose highest degree in w does not exceed 1. 

2. Unlike the _fffc-fiows, the form-invariance of z (13) under the flows generated by Hk requires extra 
restrictions on derivatives dxUi,dxUi (the string equation). 

Since 

Hk = {z'')- + l/2(^'=)o 

is a Laurent polynomial of non-positive degree, it conserves the highest degree of z{w) under the evolution 
(7), but not the lowest degree, in general. Indeed, 

Hk = z''- ((^'=)+ + l/2(^'=)o) 

is a difference between z'' and a polynomial of nonnegative degree, so that 

{Hk. z] = {-z\ z] {{{-z^)+ + l/2{z%),z} 

The second bracket in the last expression preserves the lowest degree of z. Thus the lowest degree of 
{Hk,z} equals the lowest degree of {z^^^z} = kz'^~^{z, z}. Imposing the extra restriction (10) (string 
equation), we see that the lowest degree does not exceed that of z'^~^ . Since the lowest degree of ^ is —1 
and that of 2; is —A'^, it follows that k < N + 2. 

Therefore, the form of (13) is preserved by ti, ii Toda flows with i = 1,..,N + 1, provided the string 
equation (10) holds. This completes the proof. 

It is easy to prove the converse to Proposition 1, i.e. 

Proposition 2 The string equation (10) is a necessary condition for the existence of integrable polyno- 
mial flows of the form, (9), (13), (I4) 

Proof: This essentially uses the same arguments as those leading to the result of Proposition 1. For 
the first N + 1 z flows to be consistent with the polynomial reduction it is necessary that the Laurent 

expansion of {z, z} be of the form {z, z} = f{x) + X]i>o Ui{x)/w^ . The corresponding condition for the 
z flows is {z, z] = f{x) + J2i>o Ui{x)w'\ Both can be simultaneously satisfied only if 

{z,z} = f{x) = f{x) 
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is independent of w. Differentiating with respect to U (or U), i.e. along the flow lines, using the equations 
of motion (7), and the Jacobi identity we get 

— = {Huf} = w-_ 



dti ' dx dw ' 

Since is w-dependent, while / is not, it follows that ^ = §{- — and similarly §(: = 0- Therefore 

/ is a constant in x and also constant along the flows, implying (up to a constant scaling) that the string 
equation (10) holds. 

In the following subsections we consider rational and logarithmic reductions as well as their Hamilto- 
nian structures. Polynomial solutions can then be viewed as a special case. Nevertheless the correspond- 
ing limiting procedures are rather cumbersome and it is easier to consider the polynomial case separately. 
A derivation of the canonical linearizing variables in the polynomial case is given in Appendix B. This 
closely follows the analogous procedure given in Section 4 for the logarithmic case. Local canonical 
lianearizing variables, which are constant and linear in the hierarchy times, turn out to be "Richardson's 
harmonic moments" of the Hclc-Shaw (Laplacdan growth) problem when z{w) and z(w) are identified 
with the conformal mapping from the w to the z plane. This is a simple proof of the relation (originally 
established in [12]) between the "times" of the 2dToda hierarchy and harmonic moments of the exterior 
Hclo-Shaw problem. 

An important consequence of the present section is the fact that the string equation (10) necessarily 
holds for reduced conformal mappings evolving under a sufficient number of 2dToda flows (i.e. sufficient 
for complete integrability). In the context of interface dynamics, Darcy's law thus turns out to be a 
corollary of integrability for reduced systems. 

3.2 Rational reductions. 

We now consider a rational reduction for z{w) and z{w), N > 1 given by 

zM = '-^^ = (16) 



The following Lemma states the consistency such rational reductions under the ti Toda flows for (16) 
and the U flows for (17). 

Lemma 1 The form of the function z{'w) in (16) is invariant under the ti flows and the form (17) of 
z{w~^) is invariant under the U flows for any i > 0. 

Proof. Consider the flows generated by Hk- Since 

Hk = (^'=)+ + l/2(z'=)o (18) 

is a polynomial of nonncgativc degree in w, its complement z^ — Hk is a polynomial in 1/w. The Laiirent 
expansion of brackets (8) {iJ/t, z} = —{z^ — Hk, z} around infinity therefore has the following form 

{Hk,z} = ki(T)w + ko{T) + k-i{T)w-^ + ... . (19) 

However the Lax bracket (8) also implies that it is rational of the form 

{Hk, z} = Q{w)/pN{wf. 

It follows from (19) that the degree of polynomial Q{w) does not exceed 2N + 1. 
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On the other hand, 

dt^z = P{w)/pN^{w), P{w) = PNdt^qN+l - qN+ldt^PN, 

where the degree of the polynomial P{w) also does not exceed 2N + 1. Equating the coefScients of 
polynomial P{w) — Q{w) to zero, we get a system of differential equations for r,a,b. The number of 
equations is 2N + 2. Thus wc get a compatible system of differential equations for 2N + 2 unknowns 
r, a, b. A similar argument shows the form invariance of (17) under the ik flows. 

The consistency of the rational reduction defined in (17) under the U flows and (16) under ii flows, 
respectively requires additional restrictions such as the string equation. But this only suflices to ensure 
form invariance under ti and ti flows respectively. 

Lemma 2 The string equation (10) is a sufficient condition for z{w) and z(w) ((16), (17)) to be form 
invariant under the first two-Toda flows. Flows generated by Hk, Hk with k > 1 are inconsistent with 
the rational reduction (16), (17). 

Remark: In the next section, we will show how to choose 4A'' + 2 independent commuting flows 
preserving the rational form of z (16) and z (17), generated by vector fields that are infinite linear com- 
binations of those generating the 2dToda hierarchy. 

Proof: In order to prove that z remains of the same form (16) under the ti Toda flows we would have 
to find a compatibility condition for a system of differential equations for r(T), a(T), 6(T), induced by 

(7)- 

Consider evolution of z along the flows generated by Hk- Again we write 

dz/dik = S{w)/pn'^{w), S{w) = Pjv%9jv+i - qN+idt^PN 

with S{w) a polynomial in w of order at most 2 A'' + 1. 
Since 

Hk - (z'^)- + l/2(z'^)o = hkw-'' + hk-iw-''+^ + ... + ho (20) 

is a Laurent polynomial of non-positive degree in w, the Laurent expansion of the corresponding bracket 
has the following form 

{Hk,z} = h{T)w + koiT) + k-i{T)w-^ + ... 
The definition of the Lax bracket (8) and the expansion (20) implies that 

{Hk, z} = {U{w) + R{l/w))/pN{wf 

where U{w) is a polynomial of degree at most 2N + 1 and R = {rik-i{T)w~^''^^ + . . . + rjo) is a 
polynomial in . If R{l/w) — 0, the number of equations will not exceed the number of unknowns. 
However, as we show below, the string equation does not imply the vanishing of R. 

Introduce the new variable y = w~^. The vanishing of R is equivalent to the following form of 
expansion in y 

{i?fe,z} = /o(T)+/_i(r)y-i + ... . (21) 

Consider this as a function of y and, as before, express Hk as a difference of Zk and a polynomial of 
nonegative degree. Then 

{Hk, z] = z} - {{{z'^U + l/2{z%), z] = 

= kz^-^{z, -z] - {{{z')+ + l/2{z')o), z) (22) 

Imposing the string equation (10) we see that the only flow that allows the above expansion (21) corre- 
sponds to fc = 1. A similar argument shows that only the ti flow preserves the form (17) of z even if the 
string equation is imposed. 

As seen from the above proof, there are only two flows generated by evolution operators of the form 
(9) compatible with the rational reduction. In fact, we should not expect more invariant flows associated 
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to the simple poles at w; = and w; = oo. In the polynomial case, the number of invariant flows was 
cqnal to the number of variables (polynomial coefficients), since one can associate n invariant flows to a 
pole of n th order, and the poles at zero and infinity are immovable. 

However, below we introduce additional flows related to movable singularities of z{iu) and z{w~^) 
which do preserve the rational reduction (16), (17), extending results of Krichever for the KP case. 

3.3 Additional flows for rational reductions of dKP hierarchy. 

In this subsection we recall the theory of flows related to poles at a finite number of finite points applying 
an approach previously used by Krichever for dKP hierarchy. In the dispersionless limit, the dKP and 
IdToda hierarchies are quite similar, while in 2dToda the existence of finite-dimensional reductions 
requires extra constraints. 

Let us start by recalling [5], that on the phase space of extended Benney systems, i.e. rational dKP, 
there arise some new flows related to the pole structure of the corresponding maps. These additional 
flows were introduced by Krichever (see [9]). 

Consider the partial fraction expansion under the above assumptions and the flows of the ti type 
only. For the dKP hierarchy [5], [19] these reductions have the form 



z{w) = w + Mo + y "° . (23) 

The new flows associated to the poles are deflned similarly to the polynomial case 

du,^z = {Bk,c,,z}, a = oo,l,2,... fc = 0,l,2,... (24) 
The evolution functions associated with the pole structure of z are as follows : 

Bk,oo = {ziw)'')>o, ()>o := 0+ + Go 

for an immovable pole at inflnity, while for each flnite-distance pole there appear additional flows with 
evolution functions: 

Bk^a = {ziw)'')a, -Bo, a = \og{Wa - w) . 

Here z{w)a denotes the negative part of a formal Laurent expansion of z{w) near the pole Wa- 

These additional flows commute amongst themselves and with the ordinary IdToda or dKP flows (asso- 
ciated with poles at infinity) . 

3.4 Additional invariant flows of 2dToda system. 

In the following, we will look at reductions of the 2dToda systems analogous to the above Benney- 
Krichever reductions of KP. We will thus choose the rational functions z{w), z{l/'w) appearing in (16) 
and (17) to have only simple poles and express these in the partial fraction form 

Now, introduce a new set of 4n -|- 2 evolution functions Hqj, Hij, Hi^^o, ^o,j, Sij, Hi^oo defined by 

Hk,j{w) = Bk,j{w) - ^Bk,j{w = 0), Hkjiy) = Bk,j{y) - ^Bk,j{y = 0), (27) 
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where k = 0,1, j = l..n and 

(28) 



-Bi,oo(w) = {z{w))>o, Bo,j = log{r{wj - w)) , Bi^j{w) = {z{w))j 



-81,00(2/) = {z{y))>o, Boj = log {r{wj - y)) , Bij{y) = {z{y))j, y = l/w. 

We denote the flow variables associated with Hqj as T2j, those associated with Hij as T2j-i, and 
that associated with Hi ^a as tq. The flow variables associated with Ho.j, Hij and Ho,oo are denoted 
T^2j,T2j-i,To respectively. Correspondingly, we denote the evolution functions 

ho i?i,oo = rw + uo/2, ho = Hi^ca = r/w + uo/2, 

h2j-i=H,,, = ^ + ^, -h2j-i=H,j = ^^ + ^Uj/wj, (29) 

h2j = Ho,j = log{wj -w) + l/2log{r/wj), h2j = Hqj = log{wj - l/w) + l/2log{r/wj). 

Then the following proposition holds 

Proposition 3 The An + 2 commuting Toda-Krichever flows 



dr^z = {hj,z}, df^z = {hj,z} 
dr^z = {hj,z} , df^z = {hj,z] ' 



j = 0..2n (30) 



preserve the rational form of z{w) and z{l/w) (26) (or equally (17, 16)) provided the string equation 
(10) holds. The dimension of the reduced phase space equals 4n + 3. 

We defer the proof, since this follows as a limiting case of the more general logarithmic reduction intro- 
duced in the next subsection. 

As in the polynomial case, the total number of form invariant flows preserving the string equation 
equals the dimension of the reduced phase space minus one. In what follows we show that these flows 
are Hamiltonian. Since the dimension of the phase space is odd and equals An + 3, it is, in fact a Poisson 
manifold whose symplectic leafs have dimension An + 2, which is exactly the number of commuting 
Toda-Krichever flows. 

The above result holds for a more general setting. Below we introduce a logarithmic reduction of the 
2dToda hierarchy and prove an analog of proposition 3 for logarithmic functions. Proposition 3 follows 
as a limiting case. 

3.5 Logarithmic flows. 

It is easier to prove the consistency of the rational reductions with the dynamics of the 2dToda system 
by flrst considering the more general logarithmic functions and then taking limits in which the branch 
points degenerate in pairs. Let us set 

z = r{x)w + u{x) + ^27=1 log(wj(a;) - w), 

z = r{x)w~'^ + u{x) + di \og{wi{x) - w~^), 

where ai,di are arbitrary constants, subject to the conditions 

n+l n+1 

^ai=0, ^ai = 0, (32) 

i=l i=l 

which ensure absence of logarithmic singularities at inflnity. 
Introduce evolution functions as follows: 

Ho = r{x)w + ^u(x). Ho = f{x)w~^ + ^u{x), 

Hj = \og{wj{x) - w) + ilog(r(,x)/w,(x)) j = l..n+l. (33) 

Hj = \og{wj{x) - w-^) -h \ \og{r{x)/wj{x)). 

We then have 
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Proposition 4 The string equation (10) is a necessary and sufficient condition for the existence o/2n+4 
commuting flows on the 2n + 5 dimensional space of functions of the form (31) 



d^z = {Hi, z} , df,z = {'Hi,z\ 
drtZ = {ni,z}, df,z={ni,z} 



,i = 0..n+ l...n+ 1 (34) 



In other words, the flows generated by (34), (33) are tangential to the manifold of such logarithmic 
functions if the string condition is imposed (and conversely). We therefore have 2n + 4 flows leaving 
invariant a 2n + 5 dimensional sub-manifold of the 2dToda system. These will be shown to be inflnite 
linear combinations of vector flelds generating the 2dToda flows (7). Equations (44, 45) will show how 
to express corresponding 2dToda times ti,ti in terms of the logarithmic flow parameters ri,Tj. 
Proof 

1. Commutativity : This will be shown to follow as Corollary 1 to Proposition 5 below. 

2. Consistency of (31) with equations of motion (34)'- As in the polynomial and rational cases the 
consistency of the d^z and df^z equations with the logarithmic reduction (31) follows from the fact that 
the structure of the Lax-Poisson brackets {Hi, z}, {Hi, z} is identical to the infinitesimal deformations 
in z and z induced by ones in the functions r, u, u, WijWi. However the consistency of the d^z and ^f^z 
equations with the logarithmic reduction requires an extra constraint, the string equation. First, we 
prove this for df^z. 

Differentiating z in (31) with respect to fj using the equation of motion (34), we get 

n+l o 

d^.z = wdf.r + dfU + y^^^^ = {Hi,z}. (35) 

^ W — Wj 

The left-hand side of (35) may contain any terms that are linear in w at w = oo and simple poles at 
w = Wj,j = l..n + 1. Since the Lax-Poisson bracket (8) is a bi-derivation, the singularities in w that 
may occur in {Hi, z} consist either of simple poles at w = Wj,j = l..n -|- 1, linear terms at w = oo or a 
simple pole at w = l/wj. We show that the latter is absent. 

Note that z can be represented as a sum over Hi, i = 0..n + 1 plus a w-independent function f{x); 

n + l 

z = Y,aiHi + f{x), (36) 

i=0 

where ao := 1, f{x) = \ {u{x) — Y^^^i log ^^(a;)^ . Now, using (36) we get 

Since f{x) is independent of w, the last term in (37) contains only linear terms in w plus simple poles 
at Wj,j = l..n+ 1. The term ^j^iaj{Hj,z} could, in principle, contain poles at w = 1/wj for j / i, 
but these must cancel since {Hi,z} contains no such poles. Furthermore, this term contains no pole 
at u; = 1/wi, because Hi is omitted in the sum. Since the string equation (10) holds by assumption, 
we see that the rhs of (35) only contains terms of the type induced by infinitesimal deformations in the 
functions r,u,Wj, showing that the df^z equation, together with the string equation, is compatible with 
the reduction (31). The proof for the equation is similar. 

The above proves the sufficiency of the string equation for the validity of Proposition 4. The necessity 
is proved similarly to Proposition 2. By equation (37), if the rational reduction (31) is preserved under 
the fi fiows, {z, z} can have no pole at w = 1 / Wi for any i or at zero. Similarly to preserve (31) under 
the Tj flows it can have no poles at Wi,i = l..n + 1 or oo. But since the only possible poles in {z, z} 
are at these points, we conclude that {z, z} is constant in w and hence can only be a function (say q) of 

X,T,f 

{z,z} = q 
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Differentiating the last equation with respect to n (or n) using the equation of motion (34), the Jacobi 
identity and the definition of the Lax-Poisson bracket we get 

dq dq dHi 
dTi dx dw 

Since is w-dependent, while q is not, it follows that || = ^ = (similarly ^ = 0). Therefore 

q = const and the string equation (10) holds. 

The demonstration that the reduced phase space is of dimension 2n + 5 is similar to the polynomial 
case. We view the string equation (10) together with the logarithmic reduction (31) as a set of (2n + 5) 
first order ODE's for the (2n + 5) functions r{x), u{x), u{x), Wi{x), Wi{x), i = l..n + 1 

dr ^. _ - s du _ _ . du -. _ _ . 

— = R{r,u,u,'Wi,...,wi, ...), — = U{r,u,u,wi, ...,wi, ...), — = U{r,u,u,wi, ...,wi, ...) 

dx dx dx 

— = Wi{r,u,u,wi,...,wi,...), — = Wi{r,u,u,wi,...,wi,...), i = l..n+l 

which, at least locally, determine the dependence of these functions uniquely in terms of their values 
ro, Mo, Wo, WiQ,WiQ, i = l..n + 1 at some initial value x = xq 

r = r{rQ,UQ, uq, wi^, wi,o, ■■■,x), u = u{ro,UQ, uq, i^i.o, wi^, ...,x), u = u(ro, uo, "o, w\fi--w-i_fi..., x) 

Wi = Wi{rQ, Uo, Uo, wio—, wifi..., x), Wi = iSi{ro,UQ, uq, wi^, .., iSifl—, x), i = l-n + 1 

The 2n + 5 initial values ro,uo,uo,Wio, 11)10,1 = l..n + 1 may be viewed as coordinates of the reduced 
phase space. 

Proof of Proposition 3: This follows from taking limits in which the logarithmic branch points W2i, 
W2i-i coalesce in pairs. 
Setting 

a2i-i = 1/e, a2i = -1/e W2i = t«2j-i + e^i /2g^ 
d2i-i = 1/e, a2i = -1/e W2i = W2i-i + eUj 

in (31) we get (26) in the e ^ limit. 

The evolution functions (29) , generating flows on the space of rational reduction (26), are then 
obtained as follows: 

ho = hm<;=o Ti-o, ho = hme=o T^o _ 

h2i-i = lim^^o 7(^2, - W2i-i), h2i = hm<:=o 7(^2* - ^2,-1) 

/l2i = linie=o H2i-1, /l2i-l = linie=0 '^2i-l- 

We thus obtain evolution functions for the rational reductions, where z, z have simple poles only, as a 
limiting degenerating logarithmic case. One may deduce structures related to general rational reductions 
by degenerating arbitrary numbers of logarithmic singularities at different finite points as well as at 
infinity. 



4 Poisson structure of logarithmic reductions of the 2Toda hi- 
erarchy. 

In this section we study the Poisson structure of logarithmic reductions of the 2dToda hierarchy, which 
will be shown to define finite-dimensional completely integrable systems. We find explicit expressions 
leading to a canonical Hamiltonian structure on the phase space of rational reductions of the 2dToda 
system. 

The Poisson structure of rational reductions of the IdToda hierarchy (which are infinite-dimensional, 
having no string equation constraint) is described in Appendix C. It is related to the Poisson structure 
of the Benney system considered in [14] 
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4.1 Hamiltonians, Action- Angle variables. 

Let us now introduce the following 2n + 5 functions on the reduced phase space defined by (31) with z, z 
solutions to (10), extended by the auxiliary variable x 

r,U,U,Wi, ...,Wn+l,Wi, ...iVn+l Q,Io,h, ...In+lJoJl, ...In+1, (39) 

where the new variables are defined as follows 

= -m ^^dz, Ij = ^ - , zdz, j = l..n + 1, 

Q = m T.7=o k/w, + £^ zdz. 

Using (31), we may evaluate the contour integrals in the last equation, expressing these explicitly in 
terms of the old parametrization r,u,u,Wj,Wj, j = \..n + 1 

Iq = z{l/w = 0) = w + X]"=i % In(wj), lo = z{0) = u + J^j^i o-j 

Ij = ajz{wj^) = aj (rw~^ +u + Y^^tl at ln(«)fe - wj^yj , (40) 
Ij = djz{wj'^) = dj (rwj^ +U + J^ktl "-k ln(wfe - wj^)) , j = l-n + 1, 



(41) 



Proposition 5 The functions Ii.,Ik,Q (40, 41) O'fs linearizing variables of the system (34) satisfying 

dr,Q = 0, df,Q = 0, 

drJ_k=Sjk, drJ_k=0, (42) 
drjh = 0, drjk = -6jk- 

Proof: We first calculate derivatives of each 1^ with respect to the times Tj,j = 0..n + 1. Applying 
integration by parts, we get 

2 .dik ^ ([ d i> ( \ d ([ ( '^^ '^^ d 
^ dTj dTjJ^^ dw Jw^dw\dTjJ J^^\dTjdw drj dw J 

By the equations of motion (34) and the definition of the Lax-Poisson brackets (8) this equals 

d f _dz \ ^ f dHj / dz dz dz dz\ ^ 
yji^ dw \ dTj J /u,^ dw \ dw dx dx dw J ' 

which, by the string equation (10) reduces to: 

d ( dz \ ^ f dHj ^ ^ , ( I, j = k 

^,d^Vd7j)''''-t^''' = '^' = [o, j^k. 

The first term in the Ihs vanishes, because z^ is univalent in a neighbourhood of Wfe, A; = l..n + 1 (and 
at (X) due to (32)) and the remaining integral is evaluated by substituting expression (33) for T-Lj. We 
have thus proved that 

dh ^ 
d^j='^^- 

The rest of the proposition is proved by similar computations. 
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Corollary 1 The vector fields defining (34) commute. 



Proof: Since z, z are completely determined by the new coordinates Q, I, I, infinitesimal deformations 
of the former under the Tj flows can be expressed via the chain rule through the infinitesimal deformations 
of the latter. It follows from Proposition 5 that 

_^(z\^ 
5Ti\ z dL \ z r 5fi\ z dh 



so we have 



5 5 S 5 \ f z \ _ f _9__d W z \ _ q 



Sri Stj Stj Sri J \ z J \ dli dij dij dli J \ z 



The commutativity of the other flows are seen to similarly follow from Proposition 5. 
Choosing the Poisson structure for equations (42) in a canonical way 

we may interpret I, I as canonical linearizing variables and Q as a Casimir invariant. The equations of 
motion induced by the evolution operators Hi or Hi are seen to be Hamiltonian and generated by the 
Hamiltonians /j, /j. 

The Poisson brackets {, }p in (43) are different from {, } in (8); they define a Poisson structure on 
the finite dimensional space of logarithmic reductions of the phase space satisfying the string equation, 
while the latter (the Lax-Poisson bracket) is a dispersionless limit of the commutator. 

Remark: Linearizing coordinates similar to (40) appeared in [11] in connection with the Laplacian 
Growth problem as integrals of the Laplacian Growth (string) equation, without the introduction of the 
compatible and fi flows considered here. 

The following proposition makes explicit the fact that the string equation implies the finite dimen- 
sionality of the reduced phase space. 

Proposition 6 The functions (40) form a set of integrals of the string equation (10). The string equation 
(10) implies 

^ = ^=0 ^=1 
dx dx ' dx 



Equivalently, this may be integrated to 

Q-x = Co, Ii = Ci, li = Ci, i = 0..n + 1, 

where Ci, ct are constants in x, which may be interpreted as coordinates on the 2n-|-5 dimensional reduced 
phase space. 

Proof: This is similar to the proof of proposition (5); one simply differentiates with respect to x and 
evaluates corresponding residues. 

As mentioned in the Introduction, the Laplacian growth problem is recovered by identifying z as the 
complex conjugate of z. As seen from Proposition 6, the Casimir Q, which is proportional to the area of 
D+ grows with unit speed in physical time x, while Ik, Ik are functions of the harmonic moments of the 
boundary curve. 

To obtain similar results for the rational case we take a limit as in (38). 
Corollary 2 Let z, z he of the rational reduction form (26). Then the following 2n + 3 quantities 
7o = uo - Yh^i Ui/wi, Jo = uo - X;"=i Wi/w't 

T —1,-1 \-^n Uj r 1 , I v^?* ^7 

h^-l=rw^ +Mo + Ej=i 1/^.-^, ^ h^-l=rw^ + ^ i = l..n 
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the linearizing canonical variables for the rational reduction, i.e. variables in terms of which the equations 
of motion (30) for the rational reduction (26) have the canonical form (42). 

The logarithmic or rational flows are infinite linear combinations of the 2dToda flows (7). As shown 
in Appendix B, harmonic moments of the conformal mapping z{w) are linear in the 2dToda times. 
Therefore to express ti,ti, i = 1,2.. through Ti,fi,i = 0..n+ 1 one has to evaluate the integrals 

which are functions of coordinates z,z of the reduced phase space (31), and then express z,z through 
Q,IiJi,i = Q..n+l. Then 

Mk{Q, I, I) = Mk{z{Q, I, I), z{Q, I, I)), MkiQ, I, I) = Mk{z{Q, I, I), z{Q, I, I)) (45) 

Since, modulo integration constants 

Mi = i„ M, = ti, i > 0, 

Ii=Ti, li = fi, i = 0..n+l 
and Q = X, it follows that equation (45) expresses ti, U in terms of r,, Tj. 

5 Conclusions. 

We have derived consistent finite dimensional logarithmic and rational reductions of certain flows related 
to the 2dToda hierarchy. The requirement for the consistency of these reductions was the string equation. 
Since the latter may also be viewed as the Galin-Polubarinova equation for the Hele-Shaw (Laplacian 
growth) process, this also established the latter as a compatible constraint for such reduced 2dToda flows. 

More generally, it would be of interest to determine all finite dimensional reductions of the 2dToda 
hierarchy that imply fulfilment of the Darcy law in various forms. 

Also, it would be useful to study integrable systems connected to various kinds of Hele-Shaw fiows. 
This could include boundary conditions which are more general than a point sink at infinity. For instance, 
one can consider an exterior problem with a steady and uniform viscous flow at infinity. This can describe 
the evolution of a bubble surrounded by a viscous liquid moving through a wide channel, or a bubble 
within a viscous fiow generated by a source and a sink of equal magnitude separated by a large distance. 
In such models, the constraint on z, z (5, 6) has the form 

{z{w,x),z{w~^ ,x)} = {l/w + w)r{x) 

This is an exterior analog of the interior problem with a dipole source inside the domain. In contrast to 
a monopole sink, this implies that the first harmonic moment changes with the time x, while the others 
(including area) remain unchanged. In general any combination of n-pole sources could be considered: 

{z{w,x),z{w~'^,x)} = f{x,w) 

where / is defined by the asymptotics of the hydrodynamic potential. 

6 Appendix A: String equation and Darcy law 

We recall here the derivation of the Galin-Polubarinova equation. Consider evolution of the exterior D_ 
of a simply-connected planar domain bounded by an analytic curve on ^-plane. The curve is an image of 
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the unit circle in the w plane under a conformal mapping (3). As stated in the introduction the pressure 
field P{X, Y) is constant in the interior domain and at the boundary dV^, so that at dV^ 

= dP/dx = dP/dx + VnVP = dP/dx - {VPf, 

where we used Darcy law (2) and the pressure gradient is taken in the exterior vicinity of the curve. Since 
P is harmonic in (cf (1)), it is the real part of a homomorphic function (hydrodynamic potential) 
P = Re^(^). It follows that 



= Re 



dx dz\d 



The hydrodynamic potential can be chosen as <f){z, x) = l/27rlog x) due to the logarithmic asymp- 
totic in (1). Then 

dw{z,x) dw{z,x) f dw{z,x) 



= Re\^w{z,x) ^ 

where w{z,x) stands for the map inverse to z{w,x). Using the facts that z = z{w{z,x),x); i.e., 
dz{w,x) dw{z,x) ^ dz{w,x) dw{z,x) dz{w,x) 



dw dx dx ' dz div ' 

and w = 1/w on dV^, we arrive at (4). 

7 Appendix B: Poisson structure of polynomial reductions 

Following steps similar to the proof of Proposition 5, one arrives at the following result for polynomial 
reductions (13), (14), (9) of the 2dToda system. 

Proposition 7 The following quantities 

Mk = -r^ i ZZ~''dz, Mi = -:r^ i ZZ~^dz, 



2iTTk 2iTTk Jq 

Q = ^— ( ® zdz — Q) zdz 
4t7r VJoo Jo 

are canonical linearizing variables of the (polynomially) reduced 2dToda system (13), (14), (9), (7): 

Mk{z, z) = tk + const, Mfc(z, z) = ffe + const, Q{z, z) = x + const. 
Solutions to the string equation (10) are then defined by the set of algebraic equations 

Q{z,z)=x + co, Mk{z,z) = Ck, Mk{z,z) = Ck, 
where Ck,Ck are arbitrary constants. 

Proof: Evaluating by arguments similarly to those of proposition 5, we get 

dM, __ 1 j: oHk ^ 1 i a^,-.,,. 



dtk 2iTrj dw dtj 217: j dw 

Using the fact that 

Hk = {ziwf)+ + l{z{w)>')o = ziw)" - (ziw)''). - l{z{w)'')o 
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and the negative part of the Laurent expansion does not contribute to the integral, we obtain 



dMj 1 / , ._idz{w)\ 



dtk 2iTTj /oo 9w 

Using the fact that the evolution functions Hj = {P)- + l/2{P)o generating the fj flows, are polynomials 
of nonegative degrees in 1/w, we obtain 

In the context of the Hele-Shaw problem, the first N + 1 Richardson moments [16] equal the first 

A'^ + 1 2dToda times. Indeed, in the Hele-Shaw problem z{w) has the meaning of a conformal mapping 
from w to the z plane. In the exterior problem, the exterior of the unit circle (i.e. domain \w\ > 1) is 
mapped to the exterior V- of the of the boundary curve zlw). The exterior harmonic moments are, by 
Green's theorem 

/ z~^dzdz = (p z(l/w)z(w)~^ dz. 

Since the mapping z{w) is univalent in T>^, all zeros of z{'w) are located in and we can move the 
integration contour to infinity if i > 1. Thus, we get 



— / z{l/w)z{w)~^dz = Mi 



2j7rz 

which, by proposition 7 equals U modulo an integration constant. 



8 Appendix C: Poisson structure of IdToda hierarchy. 

In this section we consider the Poisson structure of rational reductions of IdToda system. 
Recall that for the IdToda system one takes into account only the U fiows. 

dt,z = {Hi,z}, Hi = {z{wy)+ + l/2{z{wy)o, i = 0..oo (46) 

This system is bi-Hamiltonian (for general information e.g. see [3]) with two (linear and quadratic) 
compatible Poisson structures. For generic Toda system (5), the dispersionless linear Poisson brackets 
for the "field variables" Ui,i = I..00 (5) have the following form [2] 

{un{x),Um{y)}i = 2{cn + c„ - 1) [(n + m)un+m{x)S' {x - y) + mu'n_^_^{x)S{x - y)] (47) 

where 

[1 , if A;> 
Ck=\ 1/2 , if = 
[ , if < 

while the quadratic brackets are 

^ 1— n 

{unix),Um{,y)}2 = [-^{n - ■m)un{x)u'^{x) + (^{n - m + k)un+k{x)u'^_k{x) + 

fe=l 

+ku'n_^_k{x)um-k{x))] S{x -y) + [l/2(n - m)unUm + 

+ - + 2fc)u„+feU„_fe)] d{x - y). (48) 



fc-1 

As seen from the lemma 1, the rational functions 



gJV+l(w) ^ + E»=0 O'i'^^ 
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are form- invariant under all the IdToda flows (46), without any extra restriction (e.g. no string equation 
is needed). 

We obtain corresponding Poisson structures for the variables ai,bi, by using result (48), expressing 
Ui in terms of ai,bi,i = 0..N. Both the linear (47) and the quadratic (48) Poisson structures lead to 
quadratic brackets for Oj, 6, . Namely , the second Poisson structure for (49) becomes 

{ak{x), ai{y)}2 = [ ^{l + n- k)ak-n{x)ai+n{v) + nafc_„(y)a;+„(a;)) 



n=l 



+{l-N -l)ak{x)ai{y)\6'{x-y), (50) 



{hk{x),hi{.y)}2 = [^{k - I - n)bk-n{x)bi+n{y) - nbk-n{y)bi+n{x)) + 

n=l 

+ ^bk{x)bi{y)]5'{x-y), (51) 

A: — — 1 

{ak{x), bi{y)}2 = ak{x)bi{y)S'{x - y). (52) 

The first Poisson structure can be obtained from (50 - 52) with the help of a shift by a constant 

dj — >■ cii + A6i, z{w, x) z{w, x) + A 
and using a bi-Hamiltonian nature of (47), (48). 

{ak{x), ai{y)}i = [ - I - n){ak-n{x)bi+n{y) + bk-n{x)ai+n{y)) - 

n=l 

-n{ak-n{y)bi+n{x) + bk-n{y)ai+n{x))) + — Y — -h{x)ai{y) + 

+ ''^^'^^~^^ ak{xMy)]6\x-y), (53) 

{ak{x),bi{y)}i = [^{{k - I - n)bk-n{x)bi+n{y) - nbk-n{y)bi+n{x)) + 



+ ^^^ ^ bk{x)bi{y)\5'{x-y), (54) 



{hk{x)M{y)h=Q- (55) 

In all the above expressions aN+i = 1 and = if i goes beyond the range i = Q..N + 1 (and bj = if 
3 1^ 0..7V). 

These brackets form a bi-Hamiltonian structure for rational reductions of IdToda hierarchy: 

at.z = {Hi,z}i = {Hi_i,42 (56) 

with the following Hamiltonians: 

1 ^'J+i, 



Hi = / {z'^\x))odx (57) 
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